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Abstract 

Purpose The location of a phosphorus emission can strongly 
affect its expected fate in freshwater. To date, in Life Cycle 
Assessment (LCA), fate factors for phosphorus emissions 
have been derived for continents or large countries and had 
limited spatial resolution. These fate factors do not account 
sufficiently for local variations and are not applicable globally 
In this paper, fate factors for freshwater eutrophication are 
derived for phosphorus emissions to freshwater on a global 
scale with a half-degree resolution. 

Methods For this purpose, a new global fate model for 
phosphorus has been developed on a half-degree resolution. 
The removal processes taken into account are grid-specific 
advection, phosphorus retention and water use. Aggregated 
fate factors based on archetypes and on administrative units 
are presented. 

Results and discussion The derived fate factors represent the 
persistence of phosphoms in the freshwater environment. The 
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typical fate factor of phosphorus emissions to freshwater is 
10 days and can vary more than 2 orders of magnitude among 
the grid cells (the 5th and 95th percentile are 0.8 and 310 days, 
respectively). Advection is the dominant removal process of 
phosphorus in freshwater (67.5%), followed by retention 
(27.6%) and water use (4.9%). 

Conclusions The results demonstrate inclusion of information 
on the location of phosphoms emissions to freshwater can 
improve the comparative power of the fate factor implementa¬ 
tion in LCAs. The fate factors enable consistent assessment 
and comparison of freshwater eutrophication impacts at differ¬ 
ent locations across the globe. 

Keywords Differentiation • Fate factors • Freshwater 
eutrophication • Life cycle assessment • Nutrients • 
Phosphorus • Spatial 


1 Introduction 

Freshwater eutrophication is the process of increasing algal 
growth and changing species abundance and diversity in 
surface water due to an enrichment of surface water with 
nutrients (Smith 2003). Growth of phytoplankton can lead 
to reduced light penetration, increased incidence of surface 
algal scums and deoxygenation, causing fish kills and other 
effects on ecosystems and humans (Withers and Jarvie 
2008). Eutrophication has been a global concern for some 
time (OECD 1982; NRC 1992), and is likely to increase 
because of growing population and concomitant needs for 
increased food supply, agricultural land, and fertilizer use 
(Tilman et al. 2001; MEA 2005). Indeed, increased use of 
phosphoms fertilizer, relative to the pre-industrial era, has 
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already increased the global river export of phosphorus from 
8 to 22 Tg/year (Bennet et al. 2001). 

Phosphorus is often the primary limiting nutrient for 
primary production, and thus eutrophication, in freshwater 
(Schindler 1977; Carpenter et al. 1998). However, it has 
been known since the 1960s that nitrogen may also be a 
limiting nutrient (Schindler 2006). In particular, recent work 
suggests that nitrogen and even iron may be co-limiting over 
time scales relevant to biological cycles, though it is still 
likely that phosphorus may be the controlling limiting 
nutrient over scales of multiple years (Sterner 2008). 
Empirical evidence from the United States (Schindler 
2006) corroborates this generalization. Nonetheless, the rela¬ 
tionships between phosphorus supply, ambient phosphoms 
concentration in the water column and ecological response 
are complex: the movement of phosphorus through water 
bodies, related to their capacity to assimilate phosphorus, 
varies spatially and temporally (Withers and Jarvie 2008). 

Freshwater eutrophication is commonly taken into account 
as an impact category in Life Cycle Assessment (LCA). In the 
Life Cycle Impact Assessment (LCIA) phase of LCA, inven¬ 
tory emission data are translated into potential environmental 
impacts via characterization factors (Pennington et al. 2004), 
which express the fate and effect of these stressors. In the 
LCIA framework, the fate factor (FF) represents the overall 
persistence of a substance, in this case phosphoms, in the 
environment, including removal via advection, partitioning 
and water use. Ecological effect factors (EF) express the 
relationship between the ecological damage (e.g., species 
richness changes due to eutrophication) and the mass change 
of the pollutant of interest in water or another environmental 
compartment (Rosenbaum et al. 2007). 

In the past, several LCA methods have incorporated spatial 
differentiation in their fate factors (Potting and Hauschild 
2006). It has been demonstrated that exclusion of spatial 
differentiation may diminish the relevance of an impact 
assessment (Potting and Blok 1994, 1995). Conversely, 
inclusion of spatial differentiation may enhance an as¬ 
sessment by identifying locations where emissions are 
more or less impactful. For example, acidification impact 
factors in Europe may vary by 3 orders of magnitude 
(Potting et al. 1998). Eutrophication potentials were 
found to vary by up to 1.5 orders of magnitude and 
acidification potentials varied up to 3.5 orders of magni¬ 
tude in a country-level spatial model of nitrogen and 
sulfur air emissions in Europe (Huijbregts and Seppala 
2000). Finally, a spatially resolved model for emissions 
of organic chemicals in North America found a variation 
in human health impacts of up to 8 orders of magnitude 
(Humbert et al. 2009). 

For nutrients in freshwater, the methods EDIP2003 
(Potting et al. 2005) and ReCiPe (Struijs et al. 2009) 
included the nutrient fate model CARMEN (Beusen et al. 


1995), which operates at the European scale on a country 
level. In the LCA method TRACI (Norris 2003), US 
state-level differentiated fate factors were derived for N 
emissions to water using data from Fekete et al. (2000) 
and Vorosmarty et al. (2000a b), while P emissions to 
water were not spatially differentiated. The Canadian 
LUCAS method draws from EDIP2003 and TRACI for 
eutrophication, using an adapted CARMEN model for 
phosphorus fate (Toffoleto et al. 2006). To date, most 
LCA methods for freshwater eutrophication have focused 
on one continent and operate on the relatively coarse 
resolution of countries or states. However, fate factors 
for Europe and Northern America are not comparable 
and cannot be used in a multi-continent assessment. Fate 
factors for phosphorus for other continents are not avail¬ 
able, yet inventory data for a product or service may 
easily include inputs from several countries or continents, 
forcing an LCIA practioner to use characterization factors 
outside of the geographical context for which they were 
developed. 

The goal of this study is to derive spatially explicit fate 
factors for phosphorus emissions to freshwater on a world¬ 
wide scale. For the derivation of the fate factors, a new 
global fate model for phosphorus has been developed. The 
importance of the modeled phosphorus removal processes 
for the fate factors is demonstrated and the most important 
assumptions and uncertainties are discussed. 

2 Methods 

2.1 Framework 

The LCIA approach to characterizing an emission 
involves coupling an emission’s fate in the environment 
(fate factor [FF]) with the effect (effect factor [EF]) of 
that substance (see, e.g., Rosenbaum et al. 2007). The 
fate factor (FF z; , day) represents the change in mass in a 
receiving compartment cell j (AM ZJ -, kg) due to an emis¬ 
sion in compartment cell i (S h kg day^ 1 ), i.e., 

FF i j = AM iJ /S i (1) 

In the case of eutrophication, the effect factor can be 
expressed as the change in potentially disappeared fraction 
(APDF, dimensionless) of species in receiving compartment 
j over the volume (Vj, m 3 ) per unit increase in mass of 
phosphorus (Payet 2006; Struijs et al. 2011): 

11/ = ^PDF, • Vj/AMj = APDFj • Vj/ACj ■ V } 

= APBFj/ACj (2) 

Although models can have spatially differentiated effect 
factors, this discussion assumes constant effect factors, i.e., 
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EF ; =EF. It follows that the individual fate factors (FF /y , days) 
can be summed, resulting in the cumulative fate factor of an 
emission to freshwater in cell i (FF/, days), which is related to 


the fate in cell i and all receiving cells j downstream. Including 
this summation and coupling the FF and EF, the overall 
eutrophication characterization factor is expressed in Eq. 3. 
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This equation presents the general LCIA approach to char¬ 
acterization, providing context for the following discussion, 
which focuses on the development of a spatially differenti¬ 
ated fate factor for phosphorus. In the LCIA framework, if 
the effect factor is constant, a larger fate factor means a 
larger CF and thus a higher potential environmental impact. 
A larger fate factor means a higher cumulative persistence in 
the context of freshwater eutrophication, implying phospho¬ 
rus will be available longer for algae growth and the subse¬ 
quent cascade of effects caused by that growth. The 
interaction between algae blooms and cumulative persis¬ 
tence of phosphoms has not been modeled. This is in keeping 
with current LCIA practice, in which the effects captured by 
the effect factor are separate from the fate factor. Biomass 
growth is included in the removal rate of phosphorus via 
retention and settling. 

2.2 Model structure 

Fate factors of phosphorus emissions to freshwater were 
derived for a 0.5° x 0.5° grid covering the globe. At the 
global scale, advection, retention and water use were iden¬ 
tified as the most important processes affecting phosphorus 
fate in freshwater (Bennet et al. 2001; Harrison et al. 2005; 
Mayorga et al. 2010; Seitzinger et al. 2005). In each grid 
cell, each of these processes results in phosphorus removal 
and is modeled by a single rate constant, as shown by the k 
arrows in Fig. 1. The grid cells are connected through the 
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Fig. 1 Spatial relations between the model grid cells and their environ¬ 
mental compartments 


freshwater system, with advection transporting water and 
phosphorus (£ adv in Fig. 1) to downstream grid cells. Phos¬ 
phorus is also cycled between the freshwater system and the 
sediment, via the uptake and release of phosphorus by 
particles or organisms, resulting in a net removal and reten¬ 
tion (£ ret ) in the grid cells. Furthermore, both water and 
phosphorus are transferred to agricultural soil via irrigation, 
and only a share of that phosphorus is returned to the 
freshwater system by runoff. As a result, irrigation leads to 
a net removal of phosphorus from the freshwater environ¬ 
ment (£ use ). For the grid cell delineation and connectivity, 
the 0.5°x0.5° gridded river network from Vorosmarty et al. 
was used (2000a, b). The topology of their river network is 
derived from a digital elevation model and aligned with 
independent river routing data. Grid specific discharge and 
runoff data are created by reconciling empirical discharge 
data with the discharge calculated in the water balance 
model of Fekete et al. (2002). 

As shown in Eq. 3, the cumulative fate factor (FF/) for an 
emission in cell i is the sum of the fate factors for the individual 
cell of emission and of all downstream receptor grid cells 
(FF/y). The individual fate factors are the product of the per¬ 
sistence of phosphorus in freshwater in receptor cell j (r y , days) 
and of the fraction of phosphoms transported from the source 
grid cell i to receptor grid cell j (fi j, dimensionless): 

FF ; = £ FF, y = Y/UJ ■ ( 4 ) 

j j 


The persistence in grid cell j is the inverse of the sum of 
its removal rates (k adYJ -, k rQt j , k use j, day -1 ), as presented in 
Eq. 5, in which k ady j (day -1 ) represents the removal rate of 
phosphoms by the advective flow of water out of the grid 
cell, k ret j (day -1 ) represents the phosphoms retention rate, 
and k USQ j (day -1 ) is the removal rate due to water use. 



1 


^adv, j “F ^ret, /' “F k 


ret J 


use, 7 



The fraction of an emission in grid cell i that reaches 
receptor grid cell y, fj, also depends on the same removal 
rates and is the product of all fractions of phosphoms trans¬ 
ported from one grid cell to the next, from the source to 
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receptor grid cells (with / indicating the grid cell of focus 
between the source cell i to the first cell upstream of the 
receptor cell y, and f iti = 1): 


7-1 

fij =n 


k 


ad v,/ 


^adv,/ ~t~ ^-ret,/ k use,/ 



The following sections discuss the components of Eqs. 5 and 6. 


2.2.1 Advection 


The river network of Vorosmarty et al. (2000a, b) indicates 
the downstream grid cell corresponding to each cell in the 
global model. The rate of phosphorus removed by advection 
from a grid cell is equal to the rate of water removal in 
that cell. This value is the ratio between its discharge 
(Q h km 3 day -1 ) (Fekete et al. 2002) and the total water 
volume in the grid cell ( V toth km ), which is the sum of 
the volumes of lakes, reservoirs and rivers: 


k 


ad v,/ — 


Qi 

ftot ,/ 



parameters from SPARROW are derived for the US, this 
model’s calibration includes a wide diversity of streams, lake 
and climatic conditions that are also found in other continents. 
Thus, this is the most globally representative approach cur¬ 
rently available. However, care must be taken when applying 
the model to environments that are not predominant in the US, 
such as tropical ecosystems. 

7, ^ wb ’7 7, 

^retj — v ret,wb,7 

wb j 

In Eq. 8, V wbj - and V tot j (both in km ) are the total volume of 
the water body class (the subscript “wb” indicates either river, 
lake or reservoir class) and the total water volume in grid celly, 
respectively. This fraction is a volume weighting factor that 
reflects the fraction of the time the phosphoms resides in the 
water body of focus relative to its total persistence in that grid 
cell. The calculation of k ret j is explained further in Online 
Resource 1. 



2.2.3 Water use 


The river volume was calculated from the river width, 
depth and length. River width and depth were calculated by 
empirical relationships from Wollheim et al. (2006), while 
length was determined with the relation between watershed 
surface area and river length from Vorosmarty et al. 
(2000a). Details on the river volume calculation are avail¬ 
able in Online Resource 1 . Lake and reservoir volumes were 
taken from Green et al. (2004) and the International Com¬ 
mission on Large Dams (1984, 1988), respectively. For 
application of these data to fate factor calculation, the vol¬ 
ume of cells in large lakes needed to be adjusted in the 
0.5°x0.5° data set to correct for a discrepancy concerning 
low flows of edge lake cells, as the data set assigned the 
flow in such edge cells to the centerline of the water body. In 
the model presented in this manuscript, the volume of edge 
lake cells has been grouped with the volume of the center- 
line cells to ensure consistency between the flows and 
volumes of lake grid cells (Online Resource 1). 

2.2.2 Retention 

The phosphorus retention rate is governed by two main pro¬ 
cesses: the uptake of phosphoms by biomass and its adsorp¬ 
tion to suspended solids and their subsequent physical settling 
(Hejzlar et al. 2009). The overall retention rate in a grid cell 
(kretj, day -1 ) is the volume-weighted average of the removal 
rates in the separate waterbodies (k rQtwb j, day -1 ), as shown in 
Eq. 8. These removal rates were taken from the inland nutrient 
United States fate model SPARROW (Alexander et al. 2004), 
in which these rates and other input variables are calibrated to 
predict phosphoms loads in individual streams. Though 


Phosphoms is removed from the water in a grid cell when 
water is used for domestic, industrial and agricultural pur¬ 
poses. After its use, the water is eventually returned to the 
freshwater system, but the phosphoms may be retained. The 
water used for agricultural purposes was assumed to be used 
for irrigation, in which water is cycled through the soil, 
which retains a certain fraction of the phosphoms. When 
water is used for industrial and domestic purposes, a fraction 
of the phosphoms can be removed through sewage treat¬ 
ment. Global data on phosphoms removal efficiencies of 
sewage treatment plants are lacking; since applying the 
same removal efficiency globally cannot be justified, this 
process was not modeled. 

The total amount of withdrawn water can be defined 
relative to the discharge through / W ta (dimensionless) 
(Vorosmarty et al. 2000c), so that the removal rate due to 
water use k USQ J - (day -1 ) may be expressed relative to the 
removal rate by advection, k adv j (day -1 ). Because water can 
be used multiple times, ^wta can exceed 1. When the water 
use was limited by the supply, /wta,/ was set to 500, the 
maximum of the rest of the dataset. Because sewage treat¬ 
ment is not modeled, phosphoms is not removed from water 
used for domestic and industrial purposes. This is reflected 
by the term (1 - /ditwjX in which /ditw (dimensionless) is 
the share of the total water use that is used for domestic and 
industrial purposes (Vorosmarty et al. 2005). The transport 
back to the freshwater system depends on the fraction of 
total phosphoms emissions transferred from soil to fresh 
water (f so nj, dimensionless). 

^usej = /wTA, j • (1 /dITW, y) * ^advj * (1 “/soil,7) ( 9 ) 
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In Eq. 9, / SO ii is the sum of the fractions of total phosphorus 
emission to soil transferred to freshwater in the following 
forms: dissolved inorganic phosphorus (/dip), dissolved 
organic phosphorus (/dop) and particulate phosphorus 
(/pp). Tdip and f D Q p were calculated with relationships from 
Mayorga et al. (2010), which are runoff-based regressions 
for these transfer fractions. Maps of the resulting transfer 
fractions are available as Figs. 2.1 and 2.2 of Online Resource 
1. The transfer fraction of particulate phosphoms from agri¬ 
cultural soils to freshwater was not derived due to lack of data 
coverage outside of Europe. More details on the calculation 
and results of the transfer fractions are available as Online 
Resource 1. 

2.3 Aggregation of fate factors 

The cumulative fate factors can be directly applied 
when the location of the phosphorus emission is known 
at a resolution of 0.5°x0.5°. However, in applications 
such as LCA, it is currently not possible or feasible to 
specify the location of emissions at such a high resolu¬ 
tion. Therefore, fate factors should be defined for larger 
areas, such as archetypical situations or political units. 
Archetypes provide a useful distinction by which to 
categorize and group grid cells in physically meaningful 
ways, while political boundaries are often useful, e.g., 
to the LCA practitioner. 

Aggregation over larger areas can be done by emission- 
weighted averaging of the fate factors, as shown in Eq. 10: 

FF re V g e r = ^T-E FF '- £ ' 4°) 

The regional average (FF^^ e , days) is calculated 

over all grid cells with a nonzero fate factor in a region. 

is the emission of the region in grid cell / (kg day -1 ). 
Due to limited availability of freshwater phosphorus 
emission data, 0.5°x0.5° gridded population estimates 
taken from CIESIN et al. (2005) were used as a proxy 


for point emissions to freshwater. Because point emis¬ 
sions to freshwater are mainly related to human sewage 
water, population is a good proxy (Wegener Sleeswijk 
et al. 2008). 

An archetypical distinction can be made between emis¬ 
sion locations based on the presence of large water bodies 
downstream. Without large water bodies, there is often a 
relatively rapid transfer of phosphorus to the ocean. The 
choice of archetypes can be based on the cumulative distri¬ 
bution of the fate factors (Fig. 3.1 of Online Resource Al). 
This distribution is the sum of the individual distributions of 
the two archetypes, so that an inflection point in the cumu¬ 
lative distribution indicates a demarcation between one ar¬ 
chetype and another. 

3 Results 

3.1 Fate factors and aggregation 

The cumulative fate factors for phosphoms emissions to fresh¬ 
water, FF„ are shown in Fig. 2, available online in color. These 
represent the cumulative persistence of the phosphoms in the 
freshwater environment from the point of emission to the 
ocean. Note that one fifth of all grid cells have a discharge 
of zero; these are arid, and evaporation exceeds precipitation 
on a yearly basis. In the model, this implies that no water 
bodies are present in these grid cells, and FF 7 is not applicable. 
The phosphoms emitted in a non-arid grid cell that is sur¬ 
rounded by arid grid cells (such as in central Australia), is 
transported to the arid grid cell as defined by the river net¬ 
work. In those cases, FFi represents the cumulative persistence 
from its emission location to an arid grid cell. FF 7 exhibits 
considerable variability: the 5th and 95th percentiles of the FFi 
in non-arid areas are 0.8 and 310 days, respectively. The lower 
values are found in the coastal areas, where phosphoms is 
quickly transported to the oceans. The higher values are found 
in areas upstream of large lakes or reservoirs, where the 
persistence of the phosphoms is high. 


Fig. 2 Cumulative fate factors 
for phosphorus emissions to 
freshwater (FF Z , in days) on 
half-degree resolution 
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3.2 Archetypes for phosphorus persistence 

An inflection point in the ranked cumulative fate factors occurs 
around 30 days (see Fig. 3.1 of Online Resource 1); this value 
was selected to differentiate between two archetypes. The 
distinction can be interpreted as differentiating between emis¬ 
sions located upstream of large lakes and reservoirs and those 
located downstream. This is reflected in the classification of 
the fate factors in Fig. 2: darker grays (in color, blue and 
purple) reflect the emission locations upstream of large water 
bodies, with a FFi exceeding 30 days; lighter grays (in color, 
yellow and green) reflect the emission locations downstream 
of large water bodies, with a lower FF,. The FF^[^ e of these 

two archetypes are calculated according to Eq. 10, yielding 
averages of 610 days for locations upstream and 5.6 days for 
locations downstream of large water bodies. 

3.3 Regional averages for phosphorus fate 

The country averages for cumulative fate factors are presented 
in Fig. 3. A table of the country and of the continent average 
cumulative fate factors, including data on spatial variability, is 
available in Online Resource 2. Countries on the Arabian 
Peninsula are arid and thus have no fate factor. The extent of 
spatial variability is related to the size of the country and to the 
variety of hydrologic conditions in the country; small 
countries with a size of 1 grid cell exhibit no spatial variability, 
while larger countries tend toward the opposite (e.g., Russia 
has a coefficient of variation of 8.2). Out of 168 countries, two 
countries have fate factors lower than 1 day and 21 greater 
than 1 year. At the continental level, averages vary between 
4.8 days for Oceania and 474 days for Africa, with a world 
average of 130 days. 

3.4 Dominant phosphorus removal process 

The influence of the separate removal processes described in 
the Methods section was determined by calculating the net 


removal rate for each process for every grid cell. Details on the 
calculation are given in Online Resource 1. For 67.5% of the 
non-arid grid cells, the largest net removal process is 
advection. Retention is the most important process for 
27.6% of the grid cells. Most of the grid cells where 
retention is dominant are located in regions upstream of 
large water bodies with a high residence time of water, 
i.e., the North American and African Great Lakes, areas 
of Northern Canada, and around the Finnish and West 
Russian Lakes. Water use is the dominant removal pro¬ 
cess for the remainder of the grid cells (4.9%). These 
factors are found in cells located around arid areas like the 
central United States, central Australia and central Asia. The 
spatial distribution of dominant processes is shown in Fig. 4. 
The removal mechanisms and their respective consequences 
on fate factors are further detailed and discussed in Fig. 3.2 of 
Online Resource 1 . 

4 Discussion 

After comparing the current fate factor model with other 
nutrient fate models, several sources of uncertainty in the 
fate factor calculations are presented. Subsequently, possible 
advances in the impact characterization of freshwater phos¬ 
phorus are addressed. 

4.1 Comparison with other models 

This model is the first to derive fate factors for phosphorus at 
the global scale at a half-degree resolution. This development, 
however, does not facilitate comparison to existing models. 
The NEWS models (Mayorga et al. 2010), which served as a 
conceptual basis of the current model, were aimed at estimat¬ 
ing nutrient exports through rivers to the ocean, and therefore 
do not focus on the inland fate of nutrients. These models do 
take advective transport into account by employing the same 
river network; furthermore, retention and water use are both 
modeled as removal processes. SPARROW does focus on the 


Fig. 3 Population-weighted 
country averages of the 
cumulative fate factors for 
phosphorus emissions to 
freshwater 
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Fig. 4 Spatial distribution of 
dominant processes for all 
cumulative fate factors (largest 
net removal rate) 



inland nutrient fate but models freshwater nutrient yield of 
individual streams at the scale of the US, based on data 
from the US Geological Survey (Alexander et al. 2004). 
SPARROW uses a different river network and disregards 
water use. A substantial component of both models is 
dedicated to estimating emissions, while emissions would 
serve as an input in the LCIA. Therefore, fate factors 
cannot be extracted from these models. TRACI does not 
model the fate of phosphorus, assuming 100% of phos¬ 
phorus emissions will reach ecosystems sensitive to P 
(Norris 2003). 

Struijs et al. (2011) provided fate factors for phosphorus 
emissions to European freshwater by applying the nutrient 
fate model CARMEN, to which this manuscript’s model 
results for Europe can be compared. CARMEN models the 
advective transport of phosphorus on a high resolution 
(l/6°xl/6°, with 124,320 grid cells covering Europe), 
though it does not include retention and water use as re¬ 
moval processes. In Struijs et al. (2011), partial fate factors 
for phosphorus sewage (i.e., point) emissions to freshwater 
in all European countries were derived. Summation of these 
partial fate factors results in a European fate factor of 
111 days. This fate factor is approximately three times larger 
than the current model’s European average of 40 days. A 
large portion of this difference is explained by the inclu¬ 
sion of retention and water use as removal processes, 
because the fate factor increases tenfold when these 
processes are excluded. Remaining differences may be 
related to the different spatial coverage of the models or 
the averaging method. 

4.2 Uncertainties 

First, when summing the volume contributions of the three 
types of water bodies, which are relevant in the calculation 
of advection, retention and water use, the volume data of a 
water body should apply only to the grid cell of focus. For 
river and lake volumes, this is carried out. Reservoirs, 
however, often have a long shape stretched over multi¬ 
ple grid cells and multiple inflow locations in different 


grid cells, but their volume data is located in one grid 
cell (where the dam is located). This input dataset 
limitation results in an overestimation of the fate factors 
located between the beginning and the end of the res¬ 
ervoir. Any redistribution of the reservoir volume over 
upstream grid cells would introduce even more uncer¬ 
tainty and therefore has not been performed. 

Second, due to lack of global data on sewage treatment 
efficiencies, sewage treatment is not implemented in the 
calculation of the removal rate by water use, k use . In the 
future, in order to include sewage treatment in the model, 
efficiencies of sewage treatment for phosphorus could be 
derived in a manner similar to that employed by Bouwman 
et al. (2005) for nitrogen. This approach calculated the average 
sewage treatment efficiency from data on efficiencies of sev¬ 
eral sewage treatment plant types and their prevalence within 
countries. If sewage treatment were included, all fate factors 
would decrease in areas where sewage treatment is present 
and water use has a non-negligible impact on the cumulative 
fate factor. 

Third, a closed water cycle is assumed in the calculation of 
k uso . In heavily irrigated areas, however, the evaporation of 
water during the irrigation process may result in additional 
retention of phosphorus in the soil and a reduction of the 
discharge, affecting both k usc and k ady . Additional research is 
needed to incorporate these effects in a quantitative way. 

Fourth, the transfer fraction from agricultural soils to 
freshwater, relevant for k use , does not include contribution 
of particulate phosphorus to the transfer fraction. It was 
possible to calculate / soi i, including particulate phosphorus, 
for European watersheds, by using the relation between 
phosphorus emissions to soil and the transport of phospho¬ 
rus to freshwater from Beusen et al. (1995). Details of this 
calculation are given in Online Resource 1 and a map of the 
transfer fractions is shown in Fig. 4.1 of Online Resource 1. 
With the inclusion of f PP , the / soi i increase is 10% or less for 
Northern and Central Europe. For Southern Europe and the 
Norwegian coast, / soi i can be several times (typically two) 
higher due to higher erosion rates. To understand the effect 
of the potential increase of / soih FF, was computed with an 
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estimated f tot twice larger than the original / soi i. We found 
that the FF* increases by a maximum of 0.7%. It can be 
concluded that the influence of the underestimation of / soi i 
on k use is relatively small. 

4.3 Fate factors for other phosphorus emissions 

Once improved data make it possible to include particulate 
phosphorus in the transfer fraction from soil to water (/^ oil ) 
on a global scale and with a half-degree resolution, fate 
factors can be derived for phosphorus emissions to (agricul¬ 
tural) soils. In order to do so, the fate factors for emissions to 
freshwater (FF,) need to be multiplied with the improved 
f soil . Further research on improving this transfer fraction is 
important, because agricultural emissions to soil contribute 
strongly to the increased phosphorus concentrations in 
freshwater. With this addition, it will be possible to provide 
weighted averages based not only on point source emis¬ 
sions, as in this paper, but also on agriculture production 
based emissions, which have a different global distribution. 

4.4 Aggregation of fate factors 

By focusing on physically based differences in the data, the 
use of archetypes preserves the differentiation among the fate 
factors in a meaningful way. In this case, the two archetypes 
are grid cells upstream and downstream of large water bodies. 
Aggregation according to continents and countries has less 
physical meaning and introduces uncertainty through using 
population as an emission proxy, but provides LCA practi¬ 
tioners with fate factors that are easier to relate to inventory 
data. (There is thus a trade-off between differentiation and 
practicability.) Since archetypical groupings in other impact 
categories, such as freshwater ecotoxicity, may be strongly 
influenced by the cumulative residence time of water, it would 
be practical to choose the cumulative residence time of water 
instead of the phosphoms itself. Choosing identical archetyp¬ 
ical definitions would create the same regions of archetypes 
across impact categories, reducing spatial data demands for 
the LCA practitioner while providing meaningful distinctions 
between emission locations. 

4.5 Characterization factors 

The presented fate factors for phosphorus emissions to 
freshwater can be multiplied with effect factors for phos¬ 
phorus emissions to freshwater to yield characterization 
factors for freshwater eutrophication. Struijs et al. (2011) 
derived concentration dependent effect factors, with a typi¬ 
cal value of 200 PDF m 3 kg~\ This results in a 5th and 95th 
percentiles for the characterization factors of 160 and 6.2><0 4 
days PDF m 3 kg -1 , respectively. This information can be 
used in endpoint assessments to aggregate various stressors 


to an overall ecosystem quality damage score. The effect 
factor could be improved by differentiating between effects 
in rivers and lake, e.g., as done by Payet (2006). Further¬ 
more, effect factors from Payet scale with trophic state (e.g., 
oligotrophic vs. mesotrophic), implying ecosystem effects at 
all levels of P concentration. The water body specific effect 
factors could then be applied to the cumulative persistence 
of the phosphorus in the corresponding water body. Further 
research is necessary to fully distinguish lake and river 
contributions to the fate factor, and to synthesize multiple 
approaches to the effect factor. 


5 Conclusions 

The fate factors for freshwater eutrophication presented in this 
paper enable a significant increase in spatial extent and reso¬ 
lution for LCIA. Local hydrological properties have the larg¬ 
est impact on these fate factors. Characterization factors 
derived from these fate factors make it possible to consistently 
account for freshwater eutrophication impacts of globalized 
production. Previously, characterization of phosphoms emis¬ 
sions relied on factors developed in the European or North 
American context, regardless of the location of the emissions. 
A spatially explicit and consistent global approach is cmcial in 
the case of complex products whose components are manu¬ 
factured in several locations, countries, or continents. 
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